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Abstract 

The threshold electrodisintegration of the deuteron at backward angles is studied with a relativis- 
tic Hamiltonian, including a relativistic one-pion-exchange potential (OPEP) with off-shell terms 
as predicted by pseudovector coupling of pions to nucleons. The bound and scattering states are 
obtained in the center-of-mass frame, and then boosted from it to the Breit frame, where the eval- 
uation of the relevant matrix elements of the electromagnetic current operator is carried out. The 
latter includes, in addition to one-body, also two-body terms due to pion exchange, as obtained, 
consistently with the OPEP, in pseudovector pion-nucleon coupling theory. The full Lorentz struc- 
ture of these currents is retained. In order to estimate the magnitude of the relativistic effects we 
perform, for comparison, the calculation with a non-relativistic phase-equivalent Hamiltonian and 
the standard non-relativistic expressions for the one-body and two-body pion-exchange currents. 
Our results for the electrodisintegration cross section show that, in the calculations using one-body 
currents, relativistic corrections become significant (i.e., larger than 10%) only at high momentum 
transfer Q (Q 2 ~ 40fm -2 and beyond). However, the inclusion of two-body currents makes the rel- 
ativistic predictions considerably smaller than the corresponding non-relativistic results in the Q 2 
region (18-40) fm~ 2 . The calculations based on the relativistic model also confirm the inadequacy, 
already established in a non-relativistic context, of the present electromagnetic current model to 
reproduce accurately the experimental data at intermediate values of momentum transfers. 

PACS numbers: 24.10.Jv,25.10.+s,25.30.Fj 
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I. INTRODUCTION 



The problem of how to treat the relativistic dynamics of interacting, composite objects, 
such as nucleons, is highly non trivial, and a variety of different approaches have been 
developed. These fall essentially into two classes: either field-theory inspired methods, 

n n 

such as, for example, the spectator |l| and Blankenbecler-Sugar [2| covariant reductions of 
the Bethe-Salpeter equation, or methods based on relativistic Hamiltonian dynamics (for 
a review, see Ref. Q]). The former include explicitly anti-particle degrees of freedom and 
are manifestly covariant, while the latter subsume these degrees of freedom into effective 
potentials, only retain particle (i.e, positive energy) propagation in the intermediate states, 
and typically satisfy the requirements of relativistic covariance only approximately — these 
and additional issues are discussed in considerable detail in a review by Gilman and Gross {4]. 

Both these methods — field-theory inspired and relativistic Hamiltonian dynamics — have 
been used in calculations of few-nucleon properties, including binding energies, momentum 
distributions, and electromagnetic form factors. Among the many references (a large, but 
non-exhaustive, listingof them is in Ref. 4J), we only mention here the calculations of: 



deuteron form factors 



5| and triton binding energy 6| in the spectator-equation formalism; 



deuteron form factors within the framework of relativistic Hamiltonian dynamics, in the 

q n n 

front-form [7|, point-form [8J, and instant-form [9] implementations of it; binding ener gies 
and momentum distributions of A=3 and 4 nuclei in instant-form Hamiltonian dynamics 10 , 



In the present work we study the deuteron threshold electrodisintegration at backward 
angles with a relativistic Hamiltonian, including a relativistic one-pion-exchange potential 
with off-energy- shell terms as predicted by pseudovector coupling of pions to nucleons. The 
electromagnetic current is taken to consist of one- and two-body terms, the latter too derived 
from pseudovector pion-nucleon interactions. The full Lorentz structure of these currents is 
retained in the calculation of their matrix elements between the initial deuteron and final 
np continuum states. Corrections associated with the boosting of these states from the 
center-of-mass to the Breit frame, in which the evaluation of the matrix elements is actually 
carried out, are also taken into account. 

The deuteron threshold electrodisintegration proceeds predominantly via a magnetic- 
dipole transition between the bound deuteron and 1 Sq scattering state. Since the early 
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seventies, it has been known [12] that the associated (isovector) transition form factor is 
dominated, at momentum transfers of ~ (8-16) fm -2 , by the contributions of two-body 
currents of pion range. To the best of our knowledge, all calculations of the cross section for 
this process have been carried out so far within essentially a non-relativistic framework (see, 



for example, Refs. 
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dynamics study of Ref. 



4j). One exception we are aware of is the front-form Hamiltonian 



151 ]. which, however, only included single-nucleon currents. 



One of the goals of the present work is to assess the importance of relativistic effects in 
the deuteron threshold electrodisintegration. To this end, we also perform the calculation 
of the cross section with a non-relativistic Hamiltonian, phase-equivalent to the relativistic 
Hamiltonian described above, and the (standard) non-relativistic limits of the one-body and 
two-body pion-exchange current operators. 

This paper is organized into five sections. In Sec. [TT]we discuss the relativistic Hamiltonian 
from which the bound and scattering states are obtained, and the method used to boost 
these states from the center-of-mass to an arbitrary frame. In Sec. IHII we list the relativistic 
expressions adopted for the one-body and two-body pion-exchange currents, while in Sec. IIVI 
we illustrate the momentum-space evaluation of the relevant matrix elements entering the 
cross section of the deuteron threshold electrodisintegration. Finally, in Sec. |V] we present 
the results along with a discussion and concluding remarks. Details of the calculation are 
relegated in the Appendices. 



II. THE np BOUND AND SCATTERING WAVE FUNCTIONS 



The relativistic Hamiltonian used to generate the bound and scattering wave functions 



in the np rest frame is taken to be [10|, [ll|, |l6j 



H» = 2 yp 2 + m 2 + v» , (2.1) 

where consists of a short-range part vr parameterized as in the Argonne v is potential [r3] , 
and of a relativistic one-pion-exchange potential (OPEP) given by 

fnNN m fl(k) m 



E ' + Ar E 

Vi ■ p'<T 2 ■ p' <Ti -pcr 2 p x 



+ ii (E'-E) 



E' + m 



E + m 



(2.2) 
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Here m denotes the nucleoli mass, f n NN is the pion-nucleon coupling constant (f% NN /4iT 
=0.075), p and p' are the initial and final relative momenta in the center-of-mass frame, 
E = y/p' 2 + m 2 and E' = y/p' 2 + m 2 are the corresponding energies, and k = p — p' is the 
momentum transfer. The monopole form factor f n (k) = (A 2 — m 2 ) / (A 2 + k 2 ) with A n = 1.2 
GeV/c is considered in the present work. 

The //-dependent term characterizes possible off-energy-shell extensions of OPEP and, in 
particular, the value fi=l (fi = —1) is predicted by pseudovector (pseudoscalar) coupling of 



pions to nucleons, while /i=0 corresponds to the so-called "minimal non-locality" choice [18 ]. 
As shown by Friar almost three decades ago these various off-shell extensions of OPEP 
are related to each other by a unitary transformation, that is 



if terms of 2-7r-range 



H" = e -^ u H^ =0 e^ u ~ iP =0 + i/x [h^ =0 , u] , (2.3) 
(and shorter-range) are neglected. The hermitian operator U is given 



explicitly in Ref. [la]. This unitary equivalence implies that predictions for electromagnetic 
observables, such as the deuteron electrodisintegration cross section under consideration 
here, are independent of the particular off-shell extension adopted for OPEP, provided that 
the electromagnetic current operator, specifically its two-body components associated with 
pion exchange, is derived consistently with this off-shell extension. As discussed later in 
Sec. 1111} the pion-exchange two-body currents used in this work have been obtained assuming 
pseudovector coupling, and therefore the jj, — 1 prescription is taken for OPEP. From now 
on, the /i = 1 superscript is dropped from in Eq. (12. ip for simplicity. The resulting 
relativistic Hamiltonian has been constructed to be phase-equivalent to the non-relativistic 
H, based on the Argonne vi& potential. 

The momentum-space wave functions of the deuteron and np scattering states are denoted 
respectively as ^m(p; 0) and iPI_Jm s ,t{P'i 0), where p is the relative momentum and the zero 
in the argument indicates the rest frame in which the deuteron and np pair have velocit; 
V=0. The bound-state wave function with spin projection M is written as in Ref. 
whereas the wave function corresponding to a scattering state with the np pair having 
relative momentum k, and spin, spin projection, and isospin S,M$, and T (Mt=0 for np), 
respectively, is obtained from solving the Lippmann-Schwinger equation in momentum space: 

V47sM s ,t(P; 0) = 0k;5M s ,T(p; o) 



dty 

i 
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E 

ML 



dk' 1 



(2n) 3 2E k -E k 
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T^ jM ,(k,k')] k ^A4,T(p;O) , (2.4) 



where the 0's are antisymmetric two-nucleon free states — hence the factor 1/2 in the integral 
over intermediate states k' — with 



:p;0) 



y/2 



tf(k-p) 



\S+T 



5(k + p) 



5,T 
Xm s ,0 > 



(2.5) 



and E k = 2y/k 2 + m 2 and similarly for Ey. In Eq. ( 12. 4p note that the ^sa^t's satisfy 
incoming-wave boundary conditions, since these are the wave functions relevant for the 
process under consideration here, and that they have been expressed in terms of the T- 
matrix, defined as 



T b M T s>M , s (k, k') = (Vfck,T(0) I v | 0*;5Ar-r(O)) . (2.6) 

In Eq. (12.51) . Xm s ,o denotes the np spin-isospin state SMs, TMt = 0. 

Bound or scattering wave functions in a frame moving with velocity V with respect to 
the rest frame are obtained from 



^(p; V) = B(p, V) ^( P| |/ 7 , p ±; 0) = 4= 



1 - —V ■ (o-i - cr 2 ) x p 

4m 



^(pii/7>p±;0) 



(2.7) 



where 7 = 1/ \/I — V 2 , and py and pj_ denote the components of the momentum p parallel 
and perpendicular to V, respectively. Only kinematical boost corrections are retained, in 
particular the spin-dependent ones associated with Thomas precession are only included to 
order V 2 . The interaction-dependent corrections are ignored. However, it is interesting to 
note that Eqs. (12.41) and (12.71) suggest that, in order to boost the (fully interacting) scattering 
state, one only needs to know how to boost the free states. 



III. NUCLEAR ELECTROMAGNETIC CURRENT 



The electromagnetic current is taken as a sum of one- and two-body terms 

j = E Ji(Pi»Pi) +jia(p 1 / ,P2>Pi»P2) ■ (3-1) 



8=1,2 



The one-body term corresponds to the space part of the single-nucleon current jf = (J®, jj), 
with 



Jf(p/,P) 



"(Pi 



FiAQ 2 ) ^ + ^^,(Q 2 ) a^qp 



«(Pi) 



(3.2) 



where u(pi) and w(p/) (u = m^ ) are the initial and final spinors of nucleon i, a a ^ = 
(i/2) 7 Q , 7^ , and Fx^{Q 2 ) and F 2ji (Q 2 ) denote respectively the nucleon's Dirac and Pauli 



form factors, 



F a ,i{Q 2 ) = F!{Q 2 ) + F^(Q 2 ) r M /2 , a = 1, 2 . (3.3) 



These form factors are normalized as Ff(0)=F^ (0)=1 and F 2 S (0)=— 0.12 n.m. and 
F 2 (0) =3.706 n.m. (in units of nuclear magnetons). The Holder parameterization 19] of 
Fx and F 2 is used in this work. The spinor u, or rather its adjoint, is given by 

where p and E=y/p 2 + m 2 are the nucleon's momentum and energy, and Xo-t is its (two- 
component) spin-isospin state. Note that u'u—1. Finally, the four-momentum transfer g M , 
with Q 2 =~q fl q^, is taken in the Breit frame, in which the initial deuteron has momentum 
— q/2 and the final np pair has momentum +q/2, and is given by = (u,qz) with uj = 



Ef — Ei, where Ei = \jm d + q 2 /A (ma is the deuteron rest mass) and Ef = J E\ + g 2 /4 {E\. 
is the center-of-mass energy of the np pair, i.e. Ek = 2 \/k 2 + m 2 ). 

Assuming pseudovector 7r-iV coupling, the two-body current associated with pion ex- 
change is written as 

jl2(p{,P2>Pl>P2) = jl ( 2 } (Pl>P2>Pl>P2) + jl ( 2 ) (Pl>P2>Pl>P2) , (3-5) 

where j 12 is the current corresponding to the two seagull diagrams, and is the current 
associated with the pion in flight diagram. They are given by 



ji7(p'i,P2,Pi,P2) = iG v E (Q 2 )(rxxr 2 



m 2 k 2 2 — k 2 2 + m 2 
«(Pi)7T5«(Pi) [k 2 u(P2h^u(P2)] + 1 # 2 , (3.6) 

j 12 (P 1 ,P 2 ,P1,P2) - lG E {Q){T l XT 2 ) M ^ k? _ k 02 + mlk 2_ k 02 + ml 

(kx - k 2 ) [kxu(p{)j„j5u(v>x)] [^(P2)7p75«(P2)] , (3.7) 

where the four-momentum fcf = (A;?, k$), i = 1,2, has k\=E[ — Ei and k i =p' i — p«, and 
fnNN and fn(ki) are respectively the pion-nucleon coupling constant and monopole form 
factor introduced previously. The fractional momenta kx and k 2 delivered to nucleons 1 
and 2 add up to q, that is k x + k 2 =q. The nucleon isovector Sachs form factor G E {Q 2 ), 



related to F^{Q 2 ) and F 2 V {Q 2 ) by G V E (Q 2 ) = F^iQ 2 ) - (Q 2 /Am 2 )F 2 v (Q 2 ), is used in the 
two-body currents. This choice is motivated by the following considerations. In the non- 
relativistic limit, it is easy to show that the two-body pion-exchange current satisfies current 
conservation with the (non-relativistic) OPEP, obtained from Eq. (12. 2p by setting E=E'=m, 
if the same electromagnetic form factor is used in the charge operator and longitudinal 
component of the current. As shown in Appendix |A] G^(Q 2 ) is used in the non-relativistic 
expression of the charge operator. Of course, the continuity equation places no restrictions 
on the electromagnetic form factors that may be used in the transverse components of the 
current. Ignoring this ambiguity, the choice G^(Q 2 ) satisfies the "minimal" requirement for 
current conservation. In the relativistic case, we choose to keep this same electromagnetic 
form factor in the two-body currents. 

The full Lorentz structure of the one- and two-body currents is retained in the calculations 
reported here. The latter are listed, along with their respective non-relativistic limits, in 
Appendix [A] for completeness. 

Finally, in earlier published work on the form factors and threshold electrodisintegration 
cross section of the deuteron fl^.[l7] and form factors of the A=3-6 nuclei, most recently 20, 
211 ]. the contributions associated with the boosts of the initial and final wave functions were 
neglected, and only terms up to order (v /c) 2 were included in the non-relativistic expansion 
of jf, namely the well known Darwin- Foldy and spin-orbit corrections to the charge operator 
j®. Moreover, the two-body charge and current operators were taken to leading order. 

IV. CALCULATION 

In the one-photon-exchange approximation, the differential cross section for deuteron 



electrodisintegration in the laboratory frame can be expressed as 22j] 
d 2 a 



de'dfl' 



W 2 (Q 2 , q,P») + Wx{Q% q^P") tan^/2) (4.1) 



where e' and Q' are the final electron energy and solid angle, ou is the Mott cross section, and 
the invariant response functions W\ and W 2 depend on the square of the four momentum 
transfer, denoted as before by Q 2 , and the Lorentz scalar q^P^, with P^ being the four 
momentum of the deuteron in the initial state. At backward angles, the cross section above 
is dominated by Wi, i.e. transverse scattering. (Measurements of the deuteron threshold 
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electro disintegration have been performed at angles typically > 155°, see Sec. |V] for which 
tan 2 (#/2) is > 20.) Hence, in the following, we will consider only the response W\. In the 
Breit frame, defined in Sec. 1111} it can be written as 

W 1 (q,u) = ^l + q 2 /(2m d ) 2 ]T R ST (q,u) , (4.2) 

S,T=0,1 

where the contribution from the individual spin-isospin states of the final np pair is given 
by 

R ST (q,ui) = \ £ f-^-hAiqz^S^s^M^S^ + uj-Ef) . (4.3) 

6 M,Ms \ ' 

In the equation above, as in Sec. [UJ M is the spin projection of the deuteron, and S, Ms, 
and T specify the spin, spin projection, and isospin of the np scattering state, while in the 
energy-conserving ^-function and Ej are, respectively, the initial deuteron and final rap- 
pair Breit-frame energies, Ei = \jq 2 /4 + m 2 d and Ef = y 1 q 2 /4 + E 2 , with E^ = 2 a/A; 2 + m 2 . 
Finally, the three-momentum transfer q is taken along the z direction. 

The amplitude A denotes the matrix elements of the transverse components (i.e., orthog- 
onal to q) of the current operator, namely 

A(gz,k; S,M S ,T,M) = (^ Ms , T (V f )MQ^M(Vi)) . (4.4) 

Here, V , m(Vj) and iPu ; sm s ,t(^ /) are the deuteron and np scattering states boosted from 
the center-of-mass frame, where they are calculated in momentum space with the methods 
discussed in Sec. [Til to the Breit frame, in which they have velocities given by, respectively, 
V< = -(q/2)z/E i and V/ = +{q/2) Z /E f . 

By inserting Eq. ( 12.4ft into in Eq. (14.41) . the amplitude A can be conveniently decomposed 
into the sum of two terms: 

A = A PW + A FSI , (4.5) 

where 

A PW (gz,k; S,M S ,T,M) = (0 k;S M fll r(V / )|j ± ( 9 z)|^ Af (V < )) , (4.6) 

and 

A FSI (gz, k; S, M S , T, M) = £ / J^'p ' . A PW (gz, k'; S, M' s , T, M) . (4.7) 

M' s ^ ' k ~ k ' + 16 

Thus, the amplitude A PW corresponds to describing the final np states by plane waves (PW), 
while the amplitude A FS1 takes into account interaction effects in these states. 

9 



The electromagnetic current operator includes the one- and two-body terms discussed in 
the previous section. Details of the calculation of the amplitudes are reported in Appendix [El 



V. RESULTS AND CONCLUSIONS 

In this section we report the results obtained in the laboratory frame for the cross section 
of the deuteron threshold electrodisintegration at backward angles. The calculations were 
carried out with the relativistic (R) Hamiltonian of Sec. HIl including the OPEP with off- 
energy-shell extension predicted by pseudovector coupling of pions to nucleons, i.e. with 
/i=+l in Eq. (12.21) . This Hamiltonian was constructed to be phase-equivalent to the non- 



relativistic (NR) Hamiltonian, based on the Argonne vi$ potential [17j], using the methods 



111 



developed in Ref. 

In Figs. [T]and [2] we show the deuteron and np 1 S wave functions, respectively, derived 
from the R and NR Hamiltonian models (the continuum wave function is calculated at a 
center-of-mass energy of 1.5 MeV). The deuteron R D-wave is larger than the NR at inter- 
nucleon separations less than 1.5 fm, the corresponding D-state probabilities are 6.26% (R) 
and 5.76% (NR) — the difference has its origin in the local and non-local characters of the 
NR and R (/i+1) OPEP (for a discussion of this point, see Ref. [16j). However, the D- to 
S-state ratio and quadrupole moment, which are more sensitive to the wave functions in the 
asymptotic region, are respectively 0.0260 and 0.272 fm 2 in the R model, and 0.0250 and 
0.270 fm 2 in the NR [l(|. In contrast, the R and NR 1 So continuum wave functions hardly 
differ from each other, due to the vanishing of the tensor force in this channel. 

In Fig. [3] we report the results for the electrodisintegration cross section obtained in 
the laboratory frame with the R (solid line) and NR (dashed line) Hamiltonian models 
and corresponding one-body currents, given respectively in Eqs. (IA20 and flA40 . The cross 
section obtained by ignoring the boost corrections for both the initial and final states in the 
R calculation — this is equivalent to setting V=0 in Eq. (12. 7p — is displayed by the dotted 
line, labeled RNB. All calculated cross sections include the contributions of np final states 
with total angular momentum up to J=3. These contributions are responsible for filling in 
the well-known nodal structure at Q 2 ~ 12 fm -2 in the cross section obtained with one-body 



currents by retaining only the 1 So channel in the final state 12l. Il3l HQ (see Fig. [6] below). 



Finally, the inset of Fig. [3] shows the ratios of the R and RNB to the NR predictions. 
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r(fm) 



FIG. 1: The deuteron S- and D-state radial wave functions obtained with the R and NR Hamilto- 
nian models. 

The experimental data, labeled Bates, Saclay-81 and Saclay-85, are respectively from 
Refs. 



251 . |26| . |27| , and have been averaged over the interval 0-3 MeV of the recoiling np pair 
center-of-mass energy; those labeled SLAC are from Ref. |2^ |. and have been averaged over 
the interval 0-10 MeV. However, all theoretical curves in this figure, and following ones, have 
been calculated at a fixed center-of-mass energy of 1.5 MeV and at an electron scattering 
angle #=155°. The effect of the width of the energy interval above threshold of the final 
state, over which the cross section values are averaged, was studied in Ref. [l^], and found to 
be very small. The electron scattering angles in the Saclay, Bates, and SLAC measurements 
were respectively 155°, 160°, and 180°, but in fact the calculated cross section is weakly 
dependent on the specific value of the backward angle, since om tan 2 (9/2) — > ct 2 /(4e 2 ) as 
9 — > 180° (here, a is fine structure constant and e the initial electron energy). 

The inset of Fig.[3]shows that the IA (R) and IA (NR) predictions differ significantly (i.e., 
more than 10%) only for Q 2 > 45 fm~ 2 . At lower momentum transfers, the IA (R) cross 
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FIG. 2: The np So-state radial wave functions obtained with the R and NR Hamiltonian models 
at a center-of-mass energy of 1.5 MeV 

section values are within 10% of the IA (NR). Comparison of the IA (R) and IA (RNB) results 
shows the effect of the boosts corrections in the initial deuteron and final np states. We 
have verified explicitly, by switching off the Thomas precession term in the boost operator 
of Eq. (12.71) . that the dominant correction arises from the Lorentz contraction term (the 
resulting curve is essentially indistinguishable from that labeled IA (R); it is not shown to 
reduce clutter). Indeed, the IA (RNB) results can be approximately overlaid over the IA (R) 
results by multiplying the former by the factor [1 + Q 2 /(16 m 2 )}, corresponding to the square 
of the Lorentz factor 7 ~ 7^ ~ 7/ = 1/ Vl — V 2 , where V = | V| = (g/2) / ^Am 2 + g 2 /4 , and 
the deuteron binding energy and np center-of-mass energy have been neglected. A similar 
effect was discussed in Ref. [9[ in the context of a calculation of the deuteron electromagnetic 
form factors: it conforms with the naive expectation that the overlap between the initial and 
final states in configuration space is "squeezed" in the direction of motion (namely, along 
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FIG. 3: The cross sections for deuteron threshold electrodisintegration at backward angles, obtained 
with the relativistic and non-relativistic Hamiltonian models and corresponding one-body currents 



[curves labeled IA (R) and IA (NR)] and by ignoring boost corrections in the relativistic calcu 



25, 



26, 



27 



ation 



28(]. 



[curve labeled IA (RNB)], are compared with the experimental data from Refs. 
The inset displays the ratio of the IA (R) and IA (RNB) to the NR predictions. 

q) by 7 or, equivelently, that its momentum space overlap is "pushed out" by 7. 

In Fig. H] we report the cross section results obtained by including, in addition to the 
single- nucleon current, the two-body current associated with pion exchange (curves labeled 
IA+7T with NR, R, and RNB), while in the inset we display the ratios of R and RNB to NR 
predictions. For reference, we also show the IA (R) cross sections presented in Fig. [3] The 
cross section values in the R calculation are significantly smaller than those in the NR in 
the momentum transfer range Q 2 =18-40 fm -2 . In fact, close inspection of Figs. [3HH shows 
that the pion exchange contribution in the R calculation is larger than in the NR. In both of 
these calculations, this contribution is found to have the same sign, for Q 2 up to ~ 12 fm -2 , 
as the one-body contribution. At larger Q 2 values, however, the latter changes sign, and the 
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FIG. 4: Same as in Fig. [3J but with one-body and pion-exchange two-body currents. The cross 
sections obtained in the IA (R) calculation are also shown. 

resulting destructive interference between it and the two-body contribution is responsible 
for the suppression of the R cross section relative to the NR in this Q 2 region. 

To investigate the mechanisms responsible for the suppression of the R relative to NR 
predictions for Q 2 in the range 18-40 fm -2 , we have carried out two different calculations, 
the results of which are displayed in Fig. [51 In the first, labeled IA+7T (NRW), we have 
replaced the NR expressions for the one- and two-body currents with the corresponding R 
ones, in order to isolate relativistic effects in the currents. Comparison between the IA+7T 
(NRW) and IA+7T (NR) curves shows that these effects reduce the cross section, for Q 2 > 18 
fm -2 . 

In the second calculation, labeled IA+7r (NRC), we have used NR one- and two-body 
currents but R deuteron and np scattering wave functions without boost corrections — so 
this is the same as IA+7T (NR) calculation but for the replacement of the NR wave functions 
by the corresponding R ones — with the objective of isolating relativistic effects generated 
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by the Hamiltonian. As in the previous case, we find that these reduce the cross section. 

In Fig. [5] we also show the results of an R calculation in which the time components 
of the exchanged pion four-momenta in both the vertex operators and propagators of the 
two-body currents in Eqs. (IA7j) — f)A8[) are set to zero, curve labeled IA+7T (RKO). The latter 
essentially overlaps the IA+7r (R) curve. We have also verified by direct calculation that 
ignoring the retardation effects only in the pion propagators again hardly changes the IA+7T 
(R) predictions. Thus, the explicit energy dependence of the vertex operators implied by 
pseudovector coupling of pions to nucleons has a negligible effect. 




FIG. 5: (Color online) The IA+7T (R) and IA+7T (NR) predictions of Fig. H] are compared with the 
results corresponding to different approximations, labeled respectively IA+7T (RKO), IA+7T (NRW), 
and IA+7T (NRC). See text for discussion. The experimental data are from Refs. 
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28]. 



Figure M shows the results of R and NR calculations including only the 1 So channel in the 
np final state compared both to data and the R results with the "complete" np state (all 
channels up to J=3), i.e the curves labeled IA (R) and IA+7T (R) in Fig. HI The nodes at 
Q 2 w 12 fm -2 (IA) and 18 and 16 fm -2 (IA+7r) in the R and NR predictions including only 



15 




FIG. 6: (Color online)The cross sections for deuteron threshold electrodisintegration at back- 
ward angles, obtained in the IA (R) and IA+7T (R) calculations, are compared with the results of 
calculations based on the NR and R Hamiltonian models and corresponding one- and two-body 
currents, but including only the 1 Sq channel in the np final state. The experimental data are from 



Refs. 



26 



23,128] 



the 1 So channel are filled in by the contributions of higher partial waves in the complete 
calculations. The 1 So IA (R) and IA (NR) cross sections are very close to each other, and 
thus confirm the conclusions of Ref. [l5|, in which the electrodisintegration cross sections was 
calculated within a relativistic approach based on light-front-form Hamiltonian dynamics, 
including only one-body currents. 

To conclude, we find that relativistic effects in the calculations with only one-body cur- 
rents become important (larger than 10%) at momentum transfers Q 2 exceeding 40 fm -2 , 
and are due, for the most part, to boost corrections. However, when the pion-exchange 
current contributions are also taken into account, significant differences at lower Q 2 are 
obtained between the cross sections predicted within the R and NR models for the Hamilto- 
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nian and currents. The interplay between relativistic effects in the interactions and currents 
conspire to significantly reduce the cross section obtained in the R calculation in the Q 2 
range ~ 18 — 40 fm~ 2 . 

The cross section predictions based on both the R and NR models do not reproduce 
the experimental data at Q 2 > 10 fm~ 2 , thus demonstrating the inadequacy of the present 
model for the electromagnetic current operator. This conclusion corroborates that of an 
earlier (NR) study [14], and suggests the need for including additonal (short-range) two- 
body currents. 

Finally, in the present work we have not addressed the issue of current conservation within 
the R framework. Its discussion would require constructing the two-body charge operator 
associated with pion exchange, and studying the relation between the pion-exchange charge 
and current operators and the off-energy shell behavior of OPEP. This is beyond the scope 
of the present work. 
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APPENDIX A: CURRENT OPERATOR EXPRESSIONS 

In this Appendix we list the expressions for the one-body and two-body pion-exchange 
current operators. The time (charge) and space (current) components of the one-body four- 
current read, respectively: 



Pi(Pi) Pi) — N-Ni F hi + F M 



Pi Pi + icri ■ p- x pi 




+ 




2m 



( 



q ■ p, - i(Tj ■ Pi x q q ■ p. + jg-j ■ p. x q 

Ei + m E( + m 



) 



(Al) 
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Ji(pi,pi 



n'nAf u 



p, - ICTi X pi p/ + ZCTi x p/ 



Ei + m 



El + m 



2,i 



Pi - i(Ti X Pi p/ + icTi X Pi 



2m 

-i — - q x cTj 

2m 



EL- + m 



£7/ + m 



1 + 



Pi • Pi 



(El + m)(E i + m) 

Fz,i q • Pi p/ - q • p/ Pi \ 



F2,i (Ti ■ p/ q X pj + <Ti ■ pi q X p/ 



2m 



(El + m){Ei + m) 



2m (El + m)(Ei + m) 



(A2) 



where the spinor-normalization factor Ai = y (E{ + m)/(2 i%) and similarly for A/, <Ti is the 
Pauli spin operator of nucleon i, and the initial and final nucleon spin-isospin states Xo-'t! 
and XcTin are not explicitly shown. The initial and final momenta are denoted respectively 
as pi and p/, while the energy and three-momentum transfers uo and q are taken in the 
Breit frame, defined in Sec. II III after Eq. ( 13.41) . The non-relativistic limits to order (p/m) 2 



included are written as 

pf R ( P ;, Pl ) 



G 



E,i 



+ i- 



.2G 



G 



Am? 



jf R ( P ; )P 



1 + Q 2 /(4m) 2 
Ge,i , / \ . .Gma 

(Pi + Pi) + ^ — o-i X q , 



EA i 
— °i ■ Pi X Pi , 



where Q 1 



2m 2m 
and the Sachs nucleon form factors, defined as 



G 



E.i 



F\A — 



Q 2 



Am 2 

Gma = Fn + F 2i , 



2,i j 



(A3) 
(A4) 

(A5) 
(A6) 



have been introduced in Eqs. (IA3l) - flA4l) . The two-body pion-exchange currents are given 
by 

flNN fn(^) 



ji ( 2 ) (p , i,P / 2 ,Pi,P 2 ) = iKN'^NiGKQ 2 



'_kN_N 

ml 



k'i - kP + m 



Ti x r 2 



o-i + 



v q-i ■ pj ctx (o-i ■ p iy 

(J^ + m)(J5i +m) 



As 



/ °"2 • P2 . CT 2 • P2 



<x 2 ■ k 2 + 



(cr 2 ■ p' 2 ) cr 2 • k 2 (cr 2 • p 2 ) 



(E' 2 + m)(E 2 +m) 



E 2 + m E 2 + m 
+ 1 #2 , 



(A7) 



Ji? (p'i, PL Pi, Pa) = i N^N^ G V E (Q 



m 2 fcf — ki 2 + ml k 2 — k 2 2 + m 2 



ri x r 2 



(ki-k 2 ) \k\ 



°~i ■ Pi . o"i ■ Pi 



E'-y+m Ei + m 



foO I (T 2 P 2 cr 2 p 2 

2 \E 2 + m E 2 + m 



<Ti • kx + 
<x 2 • k 2 -+ 



(Ol ■ p[) O"! ■ kx (gi ■ pi) 

+m)(Ei +m) 
(<x 2 • p' 2 ) cr 2 • k 2 (cr 2 ■ p 2 



(E' 2 + m)(E 2 + m) 



•, (A8) 
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where k® = E[ — Ei and kj = p - — p«, and the product of three Pauli matrices can be further 
reduced via the identity 

Oi • Pi) (Ti (<Ti ■ Pi) = p'i(<Ti ■ pi) + Pi(<Ti ■ p-) - (T^ ■ pi) + i(p'i X Pi ) . (A9) 

To leading order in (p/m) 2 , the non-relativistic limits in Eqs. (lA7j) - (lA8l) sum up to 



j 12 (k 1; k 2 ) = -iG v E {Q 2 ) f -^f{r x xr 2 ) 



1 l CT 2 ■ k 2 j-T2 



'2 



mi 



- (k! - k 2 )( < r 1 ■ kx)(ai • k 2 ) - {^ l} 2 - 2) 2 1 . (AlO) 
APPENDIX B: CALCULATION OF AMPLITUDES 

In this Appendix we outline the method used to compute the matrix elements of the 
one-body and two-body current operators in Eqs. (I4.6p - fl4.7l) . The calculation is carried 
out in the Breit frame, in which the initial deuteron and final np pair (with center-of-mass 



energy = 2^/k 2 + m 2 ) have velocities given respectively by 

V, = q , V, = + . q . (Bl) 

2Jm 2 + 9 74 2JE 2 + q*/A 



Here the momentum transfer q is taken along the z-axis. We also define 7i=l/ y 1 — V 2 and 
similarly for jf. 

The computer codes implementing the formalism discussed below have been successfully 
tested by comparing, in a model calculation which ignored boost corrections and kept only 
the leading terms in the non-relativistic expansions of the one- and two-body currents, the 



present results with those obtained 
code. 



14| with an earlier, configuration-space version of the 



1. One-body amplitude 

Following Eq. (14. 5p . we decompose the one-body amplitude into PW and FSI amplitudes. 
The PW amplitude is written as 

Af5 body (q;k, S, M S ,T,M) = 2 J * ^^(p + q/2; V/) h,i(Pv Pi) lMP5 V.) , (B2) 
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where p is the relative momentum, p-[=3q/4 + p and pi=— q/4 + p, the factor of 2 in 
front of the integral takes into account the identical contribution coming from the current of 
nucleon 2, j_L, 2 , and finally i^m(P] Vj) and 0k;SM s ,:r(p + q/2; Vj) are the momentum-space 
deuteron and free np wave functions, boosted to the Breit frame. Making the change of 
integration variables [(py + q/2)/jf,p±] — > (p|| , P_l) , where py and p_i_ refer respectively to 

the components of p parallel and perpendicular to q, leads to the following expression for 

v 

body 



A PW 



A™ ody (q;k,S,M s ,T,M) = 2 7/ j * ^.^(p; ) B\p, V f ) j ±A (p[, Pl 



dp 

W) 

B(p, V t ) M(lfP\\ - q/2)/T« Pj (B3) 



where p{=7/p|| + q/4 + pj_ and p 1 =7/P|| — 3 q/4 + pj_, and the boost operators B(p, Vj) 
and B(p, Vj) can be read off from Eq. (12. 7p (note that under the change of variables above, 
the Thomas precession term remains unchanged, since both Vj and Vj are along q). 



It is convenient to expand the free np wave function in partial waves 23| : 



k ;5M s ,T( P ; 0) = V2 (2VT) 3 S Jt-JL £ 6 i5T [^(k)]* ^(p) xl , (B4) 

K P LJMj 

where e L5T = [1 - R L+s+r ]/2, 

ZtsMsfc) = E( LM l, SM s \JMj)Y LML (k), (B5) 

M L 

and 3^L5j are standard spin-angle functions. Inserting this expansion in Eq. ( 1B3I) gives 



Af w body (q; k, S, M s , T,M)= £ e LST Z J L f Ms (k) j[ J Mj;M (q, fc) , (B6) 

LJMj 



where 



^jM j; M(q^) =2^7/ / ^ P Xo t ^5 J J t (p)^ t (P ) V / )j^ 1 (p(,p 1 ) 

S(p,V0^[(T/P||-q/2)/T<,P±] (B7) 

and the magnitude of the relative momentum is fixed by the ^-function in Eq. ( 1B4I) to be 
|p|=fc (note that p enters in the arguments of the boost and current operators and deuteron 
wave function). For an assigned set of quantum numbers LJMj;M and ST, the function 
J(q, k) is calculated efficiently by standard Gaussian integrations over the p-directions. 
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In order to evaluate the FSI amplitude, we first introduce in Eq. (j4.7p the partial wave 
expansions for Af^ ody , Eq. (IB6I) . and for the T-matrix [23j, 

T^ M ,(k,k')=2(47r) 2 J2 ^ L '~ L e L sTe L ,sTZl^(k){z^%^ (B8) 

JMjLL' 



and then carry out the integrations over the k' solid angle to obtain 

A™ ody (q;k,S,M s ,r,M)= £ e LST Z J L f Ms (h) J S L J M/;M (q, k; FSI) 



where we have defined 

(q,fc;FSI) = £i L '- L e L ,sT 



j LJMj;M\ 



O /-oo T STJ (k k'\ 



7T JO 



-Kfc — Ek' + it 



In deriving the equations above, use has been made of the following relation: 

/ dtt k [z[f Ms (k)\ Z J L ,f Ms (k) = 5jj,5 Mj m<5ll> , 
m s J 



(B9) 



(BIO) 



(BH) 



while a standard subtraction technique [24j is employed to perform the principal value 
integration implicit in Eq. (IBlOp . 

In the partial wave expansions of the amplitudes, fully converged results for Af^body an d 
Af-iody are obtained, at the low center-of-mass energy of the final up pair of interest here 
(1.5 MeV), when all contributions with total angular momentum J < 3 are retained in the 
sum over channels. 



2. Two-body amplitude 

In this case, after rescaling the p' relative momentum as (pj|/7/,Pj_) — » (pj|>Pj_) in the 
integral, the PW amplitude reads 

Ar w bod y(q;k, s,m s ,t,m) = V J ^^4;5M S ,T(p / ;0)st( P ',v / ) 

j±,12(Pi,P2,Pl,P2) 5 (P> Vi) V>Af(P||/7i> PJ-) > ( B12 ) 

where p{=q/4 + 7/pj| +p'j_, P2= q/4 — 7/pji — Pj_ and Pj=— q/4 + p, p 2 =— q/4 — p. Rather 
than expanding the free np state in partial waves, we carry out the p' integration by inserting 
into the equation above the plane waves of Eq. (12.5j) . and obtain 

A™ body (q;k, S,M S ,T,M) = y/2 V J ^(k, V,) 

j±,12(Pl ) P 2 ' ; Pl ) P2) 5 (P ) V i)^A/(P||/7i,P±) , (B13) 
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where in the momenta p{ and p~2 the parallel and perpendicular components of the relative 
momentum p' are replaced by those corresponding to k. The three-dimensional integrations 
in Eq. (1B13[) are done by Gaussian quadratures. 

The amplitudes Af^j, od are calculated from Eq. (14.71) by direct integration over k'. To 
this end, we first reconstruct, from the channel solutions T^ J (k,k'), the full T-matrix in 
Eq. (1B8I) . by including contributions with total angular momentum up to J=3, and then 
use cubic-spline techniques to interpolate the A2^ ody , previously tabulated on a sufficiently 
coarse grid, at the k' values relevant for integration. 
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